Fermions obstruct dimensional reduction in hot QCD 
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We have studied, for the first time, screening masses obtained from glueball-like correlators in 
Quantum Chromodynamics with four light dynamical flavours of quarks in the temperature range 
1.5Tc < T < 3Tc, where Tc is the temperature at which the chiral transition occurs. We have also 
studied pion-like and sigma-like screening masses, and found that they are degenerate in the entire 
range of T. These obstruct perturbative dimensional reduction since the lowest glueball screening 
mass is heavier than them. Extrapolation of our results suggests that this obstruction may affect 
the entire range of temperature expected to be reached even at the Large Hadron Collider. 
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The quark-gluon plasma phase of QCD is currently of 
great experimental interest. The recent announcement 
that the CERN heavy ion program may have seen this 
phase, and the forthcoming start of the BNL Relativistic 
Heavy Ion Collider are just two reasons for this interest. 
The phase structure of baryon-free QCD matter is pretty 
well determined through lattice simulations |^. While 
this was sufficient in the past, experiments will soon be- 
gin to demand much more detailed information on the 
high temperature phase. In this letter we report work 
that bears on the detailed structure of this phase. 

Lattice simulations of finite temperature (T) equilib- 
rium field theories use a discretisation of the Euclidean 
formulation for partition function — 
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where is a generic field, C the Lagrangian density, and 
the Euclidean "time" runs from to l/T. The path 
integral is over Bosonic (Fermionic) field configurations 
which are periodic (anti-periodic) in Euclidean time. For 
the lattice problem, Z[T) is the trace of an appropriate 
power of the transfer matrix, T, in one of the spatial di- 
rections. All the information in the theory is encoded in 
the set of eigenvalues of T, [T) . 

Thermodynamics is determined by the largest eigen- 
value, Aq, and its derivatives with respect to various pa- 
rameters of the theory. The spectrum of screening masses 
(inverse correlation lengths), 
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is in one-to-one correspondence with the remaining eigen- 
values, and hence contain all the remaining information 
about the T > physics of the theory. Due to the lack of 
symmetry between the space and Euclidean time direc- 
tions in eq. (^), T has only a subgroup of the rotational 
symmetry of the T = Euclidean theory. Hence the spec- 
trum of /Xi is classified by the irreducible representations 



(irreps) of this new symmetry group. These symmetries 
and their irreps have been studied in [§,1). 

Dimensional reduction has provided a good qualitative 
picture of the high temperature limit of field theories. 
This reduction yields an effective theory describing equi- 
librium physics of the (3-)-l)-dimensional theory at dis- 
tance scales much greater than l/T. At these distances 
only fields with low momentum (< T) excitations sur- 
vive; field excitations with momentum scales of order T 
or more are integrated out. In practice this procedure is 
carried out perturbatively |^ . Since the Fourier modes of 
Boson fields have momenta 2mrT in the Euclidean time 
direction and Fermion modes have momenta (2n + 1)ttT 
(where n is a non- negative integer) , therefore all Fermion 
modes and all non-zero gauge field modes are integrated 
out. The result of this integration is a 3-dimensional the- 
ory of gauge fields coupled to adjoint scalars. 

All the screening masses remaining at long distances 
should then belong to the glue sector of the (3+1)- 
dimensional theory, with Fermions appearing only in loop 
corrections. This approach has had great success in deal- 
ing with the electroweak sector of the standard model 
, and has been proposed as a way of dealing with high 
temperature QCD [||,||. A key question is about the 
temperature range where this can be done reliably. 

Lattice measurements test dimensional reduction in 
three ways — 

1. The degeneracies of Hi test whether the effective 
symmetries at high temperatures are those of a the- 
ory in lower dimension. This form of dimensional 
reduction has been seen in (3-|-l)-dimensional pure 
gauge SU{2) and SU{3) theory at temperatures as 
low as 2Tr 0,|ir 



2. Perturbative matching of the couplings of the di- 
mensionally reduced theory can be tested by com- 
paring its correlators and screening masses with 
those of the full theory. Comparison of the results 
of |,|l§ with those in 3-d SU(2) Q and SU(3) 
[O gauge theories with an adjoint scalar demon- 
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strate such dimensional reduction, as do the results 
of . More recently, this reduction has also been 
tested in the gauge sector of (2 + l)-dimensional 
QCD 0. 

3. All low- lying screening masses must be shown to 
arise in the gauge sector of the theory if the per- 
turbative matching is correct. Here we report that 
(3+l)-dimensional QCD with 4 degenerate flavours 
of light dynamical quarks fails this test. 

We have performed Hybrid Monte Carlo simulations 
of SU{3) gauge theory with 4 flavours of dynamical stag- 
gered quarks at temperatures above the first-order phase 
transition temperature Tc [ p^ . The simulations were 
done at three different couplings, and the lattice size in 
the Euclidean time direction, Nr, was chosen to be 4. 
The temperature assignments at the couplings /3 = 5.1 
(T = 1.5Tc) and P = 5.15 (T = 2Tc) are made using the 
results of previous simulations at A^,- = 6 and 8 |]l6| , |l7| . 
We also made a simulation at f3 — 5.35, corresponding to 
T = (2.9 ± 0.1)Tc dTsj. We have chosen the bare quark 
mass to be m = 0.02rc- We performed a finite size scal- 
ing study at 2Tc using 4 x 16 x 10^ 4 x 24 x 10^ 4 x 24 x 12^ 
and 4 x 16^ lattices. At 1.5Tc we worked with a lattice 
size of 4 X 16 X 10^ and at 2.9Tc with 4 x 18^. The lat- 
tice sizes have been chosen so that the correlators can 
be followed in at least one direction to spatial distance 
of 2/r or more, and the transverse size of the lattice is 
large enough to avoid spatial deconfincmcnt. In thermal 
equilibrium, the HMC trajectory lengths were taken to 
be one unit of molecular dynamics time except for the 
cubic lattice at 2Tc where it was taken to be half an unit. 
Thermalisation was monitored using four variables — the 
spatial and temporal plaquettes, the Wilson line, and the 
quark condensate. 

We have analysed correlation functions constructed 
from Wilson loops of various sizes (glueball-like corre- 
lators) for the first time in full QCD at T > 0. We 
have also analysed meson- like correlators, i. e., those con- 
structed with a qq pair. The particular operators that we 
have examined belong to the scalar and vector irreps of 
the T — theory. Under the symmetry of the spatial 
transfer matrix for T > 0, i.e., the group called DJ, the 
former belong to an irrep called the Af . The vector irrep 
of T = reduces under the T > symmetry group: two 
linear combinations give the irrep and one gives an 
irrep called . 

The glueball-like operators used in this study were the 
plaquette and the planar 6-link loop. Noise reduction was 
done through a variant of the fuzzing technique. The op- 
erators were constructed at each of five levels of fuzzing, 
and a variational method was used to extract the lowest 
mass in the channel and project out the corresponding 
correlator. More details of the methodology can be found 
in H . We could follow the correlation function in the 



sector up to a distance of 4, enabling us to obtain its mass 
with a reasonable degree of confidence. 

The meson-like operators used here are local operators 
corresponding to the tt, a, p and A at T = 0. These 
screening correlators have been studied for T > in the 
A^ channel earlier [^. We report below our measure- 
ments in this channel. The first measurement in the 
channel was reported recently [|l|. At all three temper- 
atures studied here, the correlators vanished within 
errors, consistent with free field theory and unlike the 
B^ correlators in the glue sector obtained from the same 
simulations. 

We made direct measurements of the autocorrelation 
function for the 7r-like correlator, C^(2;), at T = 2Tc on 
the 4 X 24 X 10^ lattice. For the correlators at z = 
and z = 6 we found comparable autocorrelation times, 
somewhat less than 1 MD trajectory. For z > 6 autocor- 
relations are hard to measure directly. We limit them by 
the following argument. For measurements of a correla- 
tor C{z) where the relative errors e{z)/C{z) are expected 
to be similar in magnitude, strong decreases signal large 
autocorrelations. Plotting e{z)/C{z) for the 7r-like corre- 
lator as a function of z, we found no decrease for z close 
to the center of the lattice. This makes the existence 
of strong autocorrelations rather unlikely, and eliminates 
possible simulation algorithm related systematic errors 
in mass measurements. 

In the meson-like channels we have estimated the 
screening masses by two independent techniques — fit- 
ting and constructing local masses. On all lattices these 
two methods gave consistent results. 
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FIG. 1. The A'l local mass obtained using local pionic op- 
erators at T = 2Tc on four different lattices. Also shown are 
the 1-a band on the corresponding fitted masses. 

Our finite size scaling study showed that finite vol- 
ume effects are under control. Figure ^ shows the local 
masses obtained from the pion-like correlator on the four 
different lattices at 2Tc. Also shown as bands are our 
estimates of the screening masses obtained by a fit to the 
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same correlator- 



3.74 ± 0.08 (4 X 16 X 10^ lattice), 

3.92 ±0.16 (4 X 16Mattice), 

3.83 ± 0.14 (4 X 24 X 10^ lattice), 

3.76 ± 0.04 (4 X 24 X 12^ lattice). 
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Comparing the long {N^ = 24) and short {N^ = 16) lat- 
tices, we find that the effective masses stabilise very well 
to a common value. At the same time, this estimate of 
the screening mass is independent of the transverse spa- 
tial size of the lattice, indicating that the low lying eigen- 
values of the transfer matrix have become independent 
of the lattice size. Therefore we have determined the in- 
finite volume screening mass accurately. Similar lattice 
size independence of screening masses was observed in 
the other channels as well. 

The screening masses from the 7r-like and (T-like corre- 
lators were equal within errors, as were those from the p- 
like and the A-like correlators. We investigated possible 
mass splittings more carefully by making simultaneous 
fits to these pairs of correlators, taking into account the 
covariance of the data due to the fact that the correlators 
are constructed from the same configurations. This anal- 
ysis showed no hint of any mass splitting. The results 
are collected in Table |. 

In principle, the cr-like correlator in the 4-dimensional 
theory could mix with the vacuum and hence with the 
glueball-like correlator. However the 7r-like channel 
does not mix with the 4-dimensional vacuum, and the 
above degeneracy of ii.^{A^) and p.cr{A^) cannot be ex- 
plained if the latter is small due to mixing with /ig(A^). 
Hence screening masses seen in these two channels are 
characteristic of fermion-bilinear currents in the theory. 
Furthermore, our results taken in conjunction with earlier 
studies at different cutoffs indicate that vr-cr degener- 
acy and the small value of their common screening mass, 
at these temperatures, are cutoff independent. 

Our best estimates of the projection of the p-like 
screening mass is — 



Hp{A+)/T = 5.6 ±0.1 
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at T = 2Tc. For iV^ = 4, free field theory gives 
lip{Af)/T = 5.27. At other temperatures we also found 
values of iJ.p{A^)/T close to, but slightly higher than, 
this perturbative result. The ratio iig{B^)/T was found 
to be close to this, but much larger statistics are needed 
to have any reliable estimate of this number. 

In Figure |^ we show n^{A+)/T and Hg{A+)/T as a 
function of T/Tc. In contrast to the slow variation of 
fj.g{A'l)/T with changing T, we see that /i^(A^)/T is 
increasing. The latter ratio is smaller up to 2Tc, but 
there is a cross-over near 2.9Tc. It is possible that at 
sufficiently large T the perturbative value /i^(A5^)/T — 
fip{A'^)/T ~ 5.27 is reached In accordance with 

this, we tried various smooth and monotonic two param- 
eter fits to the data which have this limit as T — > cx), and 



have a finite and non-zero value at T = Tc. Two such fits 
are shown in the figure. Extrapolating these to higher 
temperatures, we find that iJ,Tr{Af) becomes compatible 
with free fermion field theory only for Tp w lOTc or more. 
It would be interesting to test this extrapolation by direct 
simulation |2^. 
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FIG. 2. The screening masses obtained using local pio- 
nic and glueball operators shown as a function of temperature 
along with phenomenological fits to the former. The fit for- 
mulae are in the same vertical order as the fit lines at large 

Our results have implications for the dimensionally re- 
duced theory. If the perturbative approach to dimen- 
sional reduction were accurate at these temperatures, 
then one should be able to integrate out the Fermions al- 
together. Their contribution would be visible only in the 
effective couplings between the zero momentum gauge 
modes. At the same time, m(^i') should be indepen- 
dent of whether it was measured in the 7r-like or p-like 
channel and, being close to the free-field value, should 
decouple from the long distance physics. However, some- 
thing non-perturbative happens in the Fermion sector at 
experimentally accessible temperatures of up to about 1 
GeV. 

Since \Xt^{A\) < iig{A^) up to 3Tc, this excitation 
cannot be integrated out of the effective theory that de- 
scribes physics at the length scale l/iJ,g{A^). This state 
of affairs cannot be accommodated into a perturbative 
description of dimensional reduction, since all Fermions 
and multi-Fermion states need to be integrated out in 
such an approach. For 3Tc < T < lOTc, /Zt(A+) still 
deviates sufficiently from its perturbative value to ob- 
struct perturbative approaches to dimensional reduction, 
although the obstruction becomes progressively smaller. 

The values oifig{A^)/T we measured in 4-fiavour QCD 
are 53% higher than in quenched QCD Q|. Although 
comparisons of glueball masses at T = in quenched 
and full QCD show much smaller shifts p3] , the increase 
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in the Debye screening mass, in going to the four-flavour 
theory is about 30% in leading order of the perturba- 
tion expansion at any T. At the same order, the in- 
crease in fig{A'^)/T is expected to be about 60%. Al- 
though perturbation theory may not be reliable here, 
this estimate shows that the large shift is not unnatu- 
ral. In contrast, ^^(A^)/T decreases by 17% when go- 
ing from quenched to 4- flavour QCD. The obstruction to 
dimensional reduction appears very clearly in 4-flavour 
QCD due to these two movements. However it is present 
even in quenched QCD, albeit in a subtle form, since 
< ^lAAt) < t^g{B+) for T « 1.5-2T,. 
This first measurement of /Ltg(A^) in full QCD thus 
sheds new light on the effective long distance theory at 
T > Tc- As far as gauge invariant correlators are con- 
cerned, the situation up to 2Tc is reminiscent less of the 
usual dimensional reduction than of the T = situ- 
ation where the long distance effective theory is built 
from scalar and pseudo-scalar fermion bilinear compos- 
ites. Thermodynamics, governed by Ao(T) alone, cannot 
be explained by this effective model. In any case, there 
are differences between the T — and T > effective 
theories. At T = chiral symmetry is broken and the 
scalar and pseudo-scalar masses are different. Here, on 
the other hand, this symmetry is restored and the masses 
are equal. 

Finally we mention that a weaker form of dimensional 
reduction may still be valid in this range of temperatures. 
This can be investigated by looking for degeneracies of all 
measurable screening masses and checking whether they 
reflect the 4-dimensional T > symmetries or an effective 
3-dimensional symmetry. This requires a measurement of 
correlations of non-local Fermion bilinear operators, and 
is left for the future. Also left for the future is a detailed 
investigation of the continuum limit of our results. A 
similar investigation in the more realistic case of 2 light 
and one medium heavy flavour of quarks is a desirable 
future computation, specially with Wilson Fermions. 
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statistics 


Afi/T 


t^n{At)/T 


t^,{At)/T 


1.5 


900 


0.04 ±0.14 


3.51 ±0.01 


3.96 ±0.04 


2.0 


1875 


0.06 ±0.24 


3.744 ± 0.001 


4.04 ± 0.03 


2.9 


825 


0.03 ± 0.06 


4.29 ±0.01 


4.16 ±0.20 



TABLE I. Statistics and measurements of various screen- 
ing masses and the splitting A/i = fj,n{A'^) — fia{Af) at dif- 
ferent T. For T = 2Tc the results come from the 4 x 24 x 10^ 
lattice. 
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